Universality classes for depinning transitions of two heteropolymers 
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We study a solid-on-solid model for depinning transitions of two directed heteropolymers with 
an effective long range interaction decaying as the inverse square of their distance. Exact calcu- 
lations of the localization length and specific heat non-universal critical exponents indicate that 
new universality classes govern the depinning transitions, different from those found for interacting 
homopolymers. Disorder driven change in the order of the phase transition and in the sign of the 
specific heat exponent are predicted. 
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The unbinding of two directed polymers 
is an important example of a class of surface depinning 
phenomena, namely surface growth j^, denaturation of 
double-stranded DNA (dsDNA) molecules Hli|, wetting 
on a substrate 0, Ii(tI , depinning of flux lines from 
a columnar defect in a type II superconductor [l^ . etc. 
Our understanding of these surface phenomena is mainly 
based on the analysis of a class of models using the solid- 
on-solid (SOS) approximation which accounts for the cru- 
cial surface fluctuations and ignores irrelevant bulk de- 
grees of freedom. A challenging topic in these systems 
is to understand whether the presence of randomness is 
relevant at the transition and thus modifies the critical 
behavior • The purpose of this paper is to discuss this 
question in a simplified model, where the original prob- 
lem of two interacting heteropolymers is formulated as 
that of an interface near a one-dimensional random sur- 
face. In the absence of long range interactions, a similar 
problem has been addressed in the past in the context of 
wetting 0, flli] and, in two dimensions, it is known 
that the transition still exists and the disorder is marginal 
in a renormalization group (RG) sense, thus requiring a 
perturbative analysis to establish whether it will move 
the system away from the pure case fixed point or not 
fill Q, 0|. A perturbative analysis to all orders 
in the strength of the disorder was carried out in 
where only a subleading order correction to the singular 
part of the free energy at the transition of the pure case 
was found. However, a perturbative expansion in an RG 
form |n| showed that disorder shifts the transition tem- 
perature and also modifies the critical behavior, though it 
was not possible to predict the new critical exponents be- 
cause the RG flux flows toward strong-disorder regions. 
Such results were confirmed by a numerical analysis in 
01, while the issue of studying the strongly disordered 
regime was addressed in ^ llfl |. However, the question 
whether the introduction of disorder gives rise to an RG 
flow toward strong-coupling regime or a new fixed point 
exists for the random system in the perturbative regime 
is still debated, as well as the problem of determining the 
order of the phase transition 
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in the presence and in the absence of long range interac- 
tions. 

We consider here an SOS model of two directed het- 
eropolymers with a long range tail interaction decaying 
as the inverse square of their distance [J 01 . This is re- 
quired to model, e.g., the effective interaction of steps on 
a vicinal surface [3, l2l| , and similar long range interac- 
tions occur in the problem of dsDNA denaturation, where 
they are meant to mimic a space-dependent stiffness, due 
to the fact that the double-stranded conformation is sig- 
nificantly more rigid than the single one 15, 19, 22, 23j. 
Modelling the original system as an interface near a wall 
and in the absence of disorder, an exact solution of the 
critical exponents P, |3| allows one to predict a non- 
trivial first order phase transition above an upper crit- 
ical dimension, second order phase transitions with non- 
universal exponents, or Kosterlitz-Thouless-like behavior 
upon changing the amplitude of the repulsive tail. When 
disorder is added along the wall, we find new universality 
classes for the depinning transitions, compute exactly the 
localization length and specific heat non-universal critical 
exponents and show that disorder can induce a change 
in the order of the phase transition as well as in the be- 
havior of the specific heat at criticality. We emphasize 
the paradigmatic nature of the present model for surface 
depinning transitions in systems where the disorder is lo- 
calized on the wall, in view of its exact solvability. 
We begin by considering a linear object directed along, 
say, the z axis but fluctuating in the transverse d- 
dimensional x space where it is subject to an isotropic 
potential V(|x|, z). The corresponding continuum Hamil- 

tonian is H = J |^ (^)^ -t- V(|a;|, 2:)|, where m is 

the line stiffness and x{z) is the relative coordinate of 
the two original directed heteropolymers. The isotropic 
potential V(|x|, z) consists of a short range pinning term 
and of a long range power law repulsive tail V(|a;|, z) — 

v°{z) Sf{x) + 1^ (27r)'*6'(|f| - A-i), where e{x) is the 
Heaviside step function, S'Kx) = A'*/^ 6{A^'^ - |f|), 
and — Sd/i'^T^Y with Sd the surface area of a d- 
dimensional unit sphere. A is a momentum cutoff set 



2 



by the short-range part of the potential. Although in 
this paper we will restrict our analysis to the case d = 1 
and s = 2, we keep d and s general in the derivation of 
our RG equations. The partition function can be evalu- 
ated by mapping the problem into a quantum-mechanics 
problem in d dimension through the transfer operator 
technique 0, 0, |2^ |2^ . In this way the most singu- 
lar part F of the line free-energy is given, in the ther- 
modynamic limit, by the lowest energy eigenvalue Eq 
of the Schrodingcr equation for the wavefunction of the 
corresponding quantum-mechanics problem. Two length 
scales (characterizing the correlations along the polymer 
and in the transverse plane, respectively) are involved in 
the depinning transition: ^|| and They are given in 
terms of the temperature T and energy Eq by the scaling 
relations ^|| = T/\Eo\ and Eq = -DT'^/mi\ (where D is 
a constant), and are thus related: — [DT/m)S,\\ [l3|. 
The phase transition can thus be studied in terms of the 
length scale alone. 

We define i;°(z) as the sum of a constant term plus 
a quenched degree of freedom along the coordinate z: 
v^{z) ~ + 5t{z). To the quenched disorder we as- 
sociate a Gaussian probability distribution P([(5r]) = 

T>[5t\ exp {^—2^ J-^ Jr^ (z)c?z^ with variance A° << 1. 
We treat the general problem of studying the whole 
set of RG equations in the presence of quenched dis- 
order by using the replica approach. After Gaussian 
integration over St, the replicated partition function is 
Zn = / X»[f°]e-i'''^ with 



T^n — 



n 



dx" 
dz 



A- 



(1) 



A^A'^d 
2Kd 



A^A'^d ^ 



where the interactions between replicas give rise to a 
renormalization of the contact energy plus an inter repli- 

cas term and we have defined Afa,b = Jq d'^x''d'^{x"' — 
x^) = or 1 for a b, depending on whether the posi- 
tions of the replicas a, b coincide inside the stripe [0, A^^] 
or not. We thus implement a momentum-shell RG ap- 
proach. In the momentum space the interaction potential 
is V„[{fcn] = ELi Va(|fc1), where 



Vad^l) = 



V^A{d,s)\k-\ 



A°A'^d 



and w = Vb - 



A"A''d 
2Ka 



2Kd 

(2) 

V, = v^diMAlll. The Integra- 
tion over fast degrees of freedom is computed perturba- 
tively up to the second order in Vq(|A:''|), using the crucial 
fact that disorder interactions are localized on the wall. 
It proves useful to introduce new dimensionless variables 



m(Z), g{l) and A(/) 

ail) - 

u{l) = B{A,d)Vo{l) 



g{l) = B{A,d)A{d,s)A'-^'^Vs{l) 
9{l) 



d — s 



A{1) = B(A,rf)^A(/) 
2Kd 



(3) 
(4) 

(5) 

and for 



where B{A,d) = ^'''^^^^ \ In the limit n 
s — d < 2, as we will assume is the case in the following, 
the small-fc behavior is governed by the following RG 
equations 



duil) 

dl 
Ml) 

dl 
dA{l) 

dl 



(2 - d)u{l) - u^{l) + g{l) - dA{l) (6) 
(2 - s)gil) (7) 

2(1 - d)A(/) - 2A(0m(0 + A2(0 (8) 



Eq.s ©-(jSJ) are to be solved subject to the initial con- 
ditions obtained by evaluating eqs.©-® at / = 0, i.e. 
with all variables set equals to their bare values. We 
stress that in the case s = 2 and d = I, the initial con- 
ditions for g and A are cutoff independent. It is easy to 
check that in the absence of disorder we find the same 
equations as in ij. Moreover, in the absence of the tail 
and in d=l, after taking the continuous limit A — > +oo 
we obtain the result that the disorder is marginally rel- 
evant as in 11]. In the remaining of this paper we shall 
consider the marginal case s = 2 only. In this case, it 
follows from eq.© that the amplitude of the long-range 
potential is scale invariant and plays the role of a pa- 
rameter in eq.©. Eq.s ©-JHIl thus have two sets of real 
solutions: the first one describes the pure system without 
disorder and 5 > Q: 



A = 



(2 - d) ± ^{d - 2)2 + 4g 



(9) 



while the second one refers to the disordered system and 
is the following: 



A* 



A* 



-d+^{d-2)^ + 4g; u* ^ (1 - d) + ^ (10) 

with g satisfying the inequality g > — — = d — 1. 

Before discussing the physics behind eqs.®-®! let us 
discuss why they can be used for an exact determination 
of critical singularities although they arise from pertur- 
bative RG. The exactness of eq.© is due to the fact that 
the thin-shell integration is analytic and therefore does 
not renormalize a singular Vg-type interaction, whose RG 
equation thus follows from dimensional analysis. Regard- 
ing eqs.® and we will now take advantage from 
the equivalence between the present system and the one- 
particle quantum-mechanics problem in d dimensions. 
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and prove that those equations correspond to the rig- 
orous solution of the Schrodinger equation at large dis- 
tances for the wave function of the state with zero energy. 
In fact, considering a radially symmetric wave function 
and defined R{r) as its radial part, the state with zero 
energy obeys the following equation 

R"+'^^R'-'^R^O (11) 

We then make the following ansatz for the function R{r) 

f u(r) - A(r) , , , 

R = const exp / ^ ^—^dr (12) 

Then, introducing (|12|) into ec ljllll and using © we 
obtain for u(r) the same equation as in ®, provided 
/ = ln{r/a) where a is a cut-off defined by the equiv- 
alence 2mVs/T^a'-^ = {2mKaA/T^A^^'')Vs. All the 
informations regarding the issue of matching the long- 
range solution of the wave function with its short-range 
part are encoded into the initial conditions. This proves 
that eqs.jni,® give us the spatial behavior of the wave 
function with E = 0. The number of zeros of R{r) co- 
incides with the number of lower {E < 0) discrete levels 
{"oscillation theorem"). In the proximity of the phase 
transition the discrete spectrum disappears and the po- 
sition of the zero of the E = wave function occurs at 
larger and larger distances from the origin. The scaling 
argument proposed in Q allows one to identify the po- 
sition of this zero with the scale £,±. Hence we need to 
find a finite scale I* at which the R{r) vanishes and then 
study how the corresponding localization length diverges 
on approaching the phase transition. 
In eq. © , the symbols s and u are used for stable and un- 
stable fixed points, respectively. Us describes the prop- 
erties of the unbound phase, while u„ is the unstable 
fixed point where a binding-unbinding transition occurs. 
For {d — 2)^ -I- Ag < d only those solutions exist while 
the one in ifTUI) becomes unphysical, because the pa- 
rameter A can be only positive. In this case the de- 
pinning transition occurs at Uu and there is a scale 
for which u{l l^) ~oo, indicating the formation 
of a bound state. To the corresponding diverging spa- 
tial scale = A^^e' it can be associated the expo- 
nent i'± — {us — Uu)^^- Then, using the scaling relation 
a = 2 — 2di>±, the specific heat exponent is also ob- 
tained. One finds that for Us — Uu > 2 there is a region 
of the phase space where the transition is first order with 
a jump of the first derivative of the free energy, while 
for < Us — u„ < 2 the transition is second order and 
the critical exponents depend upon g and are thus non- 
universal 

A stability analysis of the RG equation ©-© shows 
however that the "pure" fixed point (A = , u = u„) is 
unstable after introducing disorder. In the presence of 
quenched disorder, depending on the initial conditions. 



RG flux flows either toward strong coupling regime or 
toward the new "random" fixed point given by ea. H10() 
which is the new locus of the depinning transition, pro- 
vided (d - 1) < g < {d + 2) (3d - 2)/4. A cartoon of 
the RG flow is shown in FigQ] for d = 1 and g = 0.5: 
the presence of a new unstable fixed point for the disor- 
dered system together with the two critical points of the 
pure system is apparent. In terms of what happens into 
the argument of the function R{r) defined in ea. (|12|l . the 
presence of a "random" fixed point corresponds to the 
case u{r) — > —oo and A(r) = A"^. The phase transition 
occurring at the "random" fixed point belongs to a new 
universality class whose exponents are non universal and 
depend upon the value A*. It is possible to compute 
the exponents ly^ and a by fixing the value of A = A"^ 
in eq.lO and proceeding as in the pure problem. One 
finds that such exponents are given by the same formu- 
las as in the pure case but with a renormalized parameter 
g^i = g — dA* . A rich scenario arises where the order of 
the phase transition is changed by disorder. The range 
of parameters where a continuous phase transition occurs 
in the disordered system overlaps the one of continuous 
phase transitions in the pure system in c? = 1, but in 
higher dimensions regions of the phase space exist where 
a continuous phase transition substitutes an otherwise 
first order phase transition, thus shifting the upper crit- 
ical dimensionality toward higher values of g. 
Another important effect of disorder is to change the sign 
of the specific heat exponent. In c? = 1 and <? > 0, the 
"random" critical exponent of the specific heat is neg- 
ative, unlike that of the pure system: the specific heat 
is thus finite at the critical point unlike the divergent 
specific heat of the corresponding pure case. Similar fea- 
tures are predicted in higher dimensions where a finite 
specific heat takes the place of an otherwise divergent 
one in the range of 5- values where the "random" solu- 
tion exists. The critical exponents and a both at the 
pure and at the "random" critical point are collected in 
Table 1. 

TABLE I: Critical exponents at the pure (1) and random crit- 
ical point (2), with Qr ~ g — dA* . 

Exponents Formulas 

y(d-2)2+45 
2[l-d((d-2)2+4g)-5] 
^lU ^{d-2y+4gR 
a(2) 2[l-d((d-2)2+4gfl)~*] 



We conclude discussing the case of g = and d = 1, 
corresponding to the situation where no long range in- 
teraction is considered. The RG flow is shown in FigI21 
it is evident that the location of the "random" critical 
point coincides with that of the pure case, i.e. A* = 
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FIG. 1: A cartoon of the RG flow in d = 1 and g = 0.5. The 
disorder variance is along the y-axis and u is along the x-axis. 
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FIG. 2: A cartoon of the RG flow in d = 1 and g = 0. The 
disorder variance is along the y-axis and u is along the x-axis. 



u* = 0. Thca'C' arc, however, interesting features of the 
phase transition in the disordered system to be stressed. 
Although in the continuous Umit (A +oo) u{Q) = 
and thus the introduction of small disorder gives rise 
to an RG flow toward strong coupling regimes, at fi- 
nite cutoffs and depending on the initial condition A*^ 
and u{0), the RG flux can flow either to the strong cou- 
pling regime or toward the "p7ire" stable fixed point. For 
< u(0) = A° << 1 a separatrix exits between these 
two regions and along it the RG flux flows toward the 
critical point A* = u* = 0. 

To summarize, this paper contains an RG analysis of de- 
pinning transitions of two directed heteropolymers in the 
presence of an effective long range interaction. The ex- 
istence of a new "random" critical point has been pre- 
dicted: two exponents have been exactly computed and 
found to depend upon disorder. Depending on the ini- 
tial value of the disorder variance A", the RG flux either 
flows toward strong coupling regime or to the new "ran- 
dom" fixed point. The order of the phase transition is 
changed by disorder as well as the sign of the specific 
heat exponent. Similar results are expected to hold in 
related problems. 
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